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D . Abstract. We investigate the zero-temperature phase diagrams of the bilayer square-lattice 

Kugel-Khomskii (d 9 ) model involving entangled and singlet phases using mean- field cluster 
approach. This diagram includes interlayer singlet phase observed in fluoride K3CU2F7 and 
exotic entangled spin-orbital phases. For a monolayer case, realized in K2CUF4, we perform 
similar calculations in finite temperature and show that the alternating-orbital ferromagnet 
decays first to an entangled uniform ferromagnet and then to a paramagnetic phase. 
Published in: J. Phys.: Con]. Series 391, 012085 (2012). 
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It has been shown that quantum fluctuations are enhanced near the orbital degeneracy and 
could suppress long-range order in the Kugel-Khomskii (KK) model [lj, called below the d 9 
model. This model was introduced long ago for a perovskite KCUF3 [2], a strongly correlated 
system with a single hole within degenerate e g orbitals at each Cu 2+ ion. Kugel and Khomskii 
showed that orbital order can be stabilized by a purely electronic superexchange mechanism. 
This happens for strongly frustrated orbital superexchange |3j, and columnar Ising-type of order 
is obtained [4| in the two-dimensional quantum compass model. This model exhibits nontrivial 
£S| I symmetry properties which may be employed to perform efficient calculations for square compass 

clusters [5]. 

Orbital order occurs in a number of compounds with active orbital degrees of freedom, 
where strong Coulomb interaction localizes electrons (or holes) and gives rise to spin-orbital 
superexchange [6]. When spin and orbital pseudospins couple to each other, their order is 
^ \ usually complementary — alternating orbital (AO) order accompanies ferromagnetic (FM) spin 

order, and ferro-orbital (FO) order coexists with antiferromagetic (AF) spin order. However, 
the above Goodenough-Kanamori rules, see also [3], are not satisfied in cases when spin-orbital 
entanglement (SOE) dominates [7J, as in the spin-orbital d 1 model on a triangular lattice [8]. 

The spin-orbital superexchange KK model for Cu 2+ (d 9 ) ions in KCUF3 with 5 = 1/2 spins 
and e g orbitals described by r = 1/2 pseudospin was derived from the degenerate Hubbard 
Hamiltonian with hopping t, intraorbital Coulomb interaction U and Hund's exchange Jh 
[9]. It describes the Heisenberg SU(2) spin interactions coupled to the orbital problem by 
superexchange J = At 2 /U, 



« = "^E { n W> + ^ n fe>) (i " tM) + ^ Q - rfj (I - if) } - E z £ rf . 
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where {ri,V2, r%} depend on rj = Jh/U [9], and 7 = a, 6, c is the bond direction. In a bilayer two 
ab planes are connected by interlayer bonds along the c axis |10j (a monolayer has only bonds 
within a single ab plane). Here 

\ 3 

= 2 ~ s« • s«+ 7 , = - + Sj • s i+7 , (2) 

are projection operators on a triplet (singlet) configuration on a bond (ij), and are the orbital 
operators for bond direction 7 = a, 6, a They are defined in terms of Pauli matrices {erf , af } as 
follows: 

^^HtfiVSof), *f = iof. (3) 

Finally, J5 a is the crystal-field splitting which favors either x = x 2 — y 2 (if E z > 0) or z = 3z 2 — r 2 
(if E z < 0) orbitals occupied by holes. Thus the model Eq. (1) depends on two parameters: 
E z /J and 77. 

The spin-orbital model Eq. (JT]) describes also CuC>2 planes in La2Cu04, where indeed U 3> t 
and large E z /Jh ~ 0.27 favors holes within x orbitals [9]. The superexchange between Cu 2+ 
ions ~ 0.127 eV reproduces there the experimental value. In this paper we consider the model 
Eq. flU for K3CU2F7 bilayer and K2CUF4 monolayer compound where nearly degenerate e g 
orbitals are expected. It has been shown that the magnetic state of K3CU2F7 is described by 
the interlayer valence bond (VB) phase [11] with AO configuration whereas K2CUF4 undergoes 
a pressure driven phase transition from the FM phase with alternating orbitals at low pressure 
to the AF phase with with x orbitals uniformly occupied. 

The simplest approach is a single-site mean field (MF) approximation applied to the model 
Eq. (P). It excludes any spin fluctuations so the spin projectors n*^j can ^ e replaced 

by their mean values depending on the assumed magnetic order. In the orbital sector we apply 
then the MF decoupling for the products {t/t^L,} along the axis 7: 

t7t7 ±7 ~ (t?)^ + t?(t^) - (t7)(t7 ±7 ). (4) 

As order parameters we take t a = (rf ) and t c = (rf ) for a chosen site i = 1 (which is sufficient in 
orbital sector as t b = —t a — t c ) and we assume two orbital sublattices: each neighbor of the site 

% is rotated by tt/2 in the ab plane meaning that (r^j) = t b ( a \ The self-consistency equations 
can be solved analytically and the phase diagram can be obtained by comparing the ground 
state energies for different points in the (E z /J,rj) plane (see Ref. [10j). One finds two classes of 
solutions: (i) uniform orbital configurations (t c = ±1/2, t a ^ = =Fl/4) for global FO order, and 
(ii) nontrivial AO order with orbitals staggering from site to site in ab planes. 

In a better cluster MF (or Bethe-Peierls- Weiss) approach, introduced to capture the effects 
of quantum fluctuations, one divides the bilayer square lattice into separate cubes containing 8 
sites each and treats the bonds inside a cube exactly, and the bonds connecting different cubes in 
MF. This approach has at least three advantages over the single-site MF: (i) spins can fluctuate, 
(ii) elementary cell can double, and (Hi) we can have independent spin-orbital order parameter. 
The MF leads in a cluster to five order parameters: magnetic s = (Sf), orbital t a ^ = (r"' ), 
and spin-orbital v a ^ = (Sfri' ). The self-consistency equations can be solved only numerically 
by iterative Lanczos diagonalization of a cluster followed by the update of the mean fields (for 
details see Ref. [10]). To capture the effect of the SOE we first obtain the phase diagram 
assuming v y ' = st a ' b (factorizable SO mean field)- see Fig. H^a) and then do the calculation 
with true v a ^ to get Fig. H{b). 

Including spin fluctuations in this approach stabilizes three valence bond phases absent in the 
single-site approach. These are VBz, VBm and PVB phases, see Fig. HJa). The VBz replaces 
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Figure 1. Cluster MF phase diagrams of the d 9 bilayer model (1) obtained with: (a) 
factorizable SO mean field, and (b) independent, potentially non-factorizable SO mean field. 
Light shaded (yellow) areas mark singlet phases with spin disorder and dark (orange) shading 
indicates phases with SO entanglement. 



ordinary antiferromagnet in the negative E z region of the phase diagram and involve interlayer 
singlets accompanied by FOz configuration which can smoothly evolve for growing rj towards 
AO configuration in the VBm being the phase observed in K3CU2F7 by Manaka et al. Here 
we explain it for realistic rj ~ 0.14. In the PVB (plaquette VB) phase spin singlets are pointing 
uniformly in a or b direction within the cluster and the clusters form a checkerboard pattern so 
the unit cell is doubled. Other phases, i.e. FM, G-AFx and ^4-AF, exhibit long range magnetic 
order and can be obtained in ordinary MF approach (see Rcf. |10j). 

A different class of phases stems from the non-factorizable SO mean field and involves SOE 
- these are the ESO, EPVB and PVB-AF phases shown in Fig. QJb) . The PVB-AF phase 
connects PVB and G-AF phases by second order phase transitions and is characterized by 
fast changes in orbital order and appearance of global magnetization. The ESO phase has no 
magnetization and weak FO order. When E z grows, the ESO phase does change continuously 
into the EPVB configuration, being an entangled precursor of the PVB phase with highly non- 
uniform AF magnetization. As shown in Ref. |10| . both ESO and EPVB are characterized by 
strong on-site SOE defined as r a ' b = v a,b — st a ' b . 

Using the same cluster MF approach as above one can easily study the properties of the KK 
model for a single layer at zero and finite temperature T. At T = one finds the phase diagram 
of the form shown in the Fig. [2Ja). In the low- 77 region it involves only AF and PVB phases 
and the VBz is replaced for AFz in the absence of the second layer. The PVB area is strongly 
elongated in the horizontal directions and caps AFz down to E z < —2J where surprizingly one 
finds a narrow stripe of a paramagnetic (PM) phase. This stripe separates PVB from FMz 
phase (FM phase with FOz order), and shows that frustration in a single layer system is bigger 
than in a bilayer one where paramagnetizm is absent in T = 0. 

In the Fig. E)^b) we show the thermal evolution of spin and orbital order parameters and 
SOE for a point (E z = — 2.1J, r\ = 0.19) lying in the FM area of the monolayer phase diagram. 
This evolution involves passage through the FMz phase and ends in the PM phase. The FM 
magnetization s exhibits critical behavior staying close to 1/2 in the FM phase and dropping 
quickly to zero at the FMz-PM phase transition. The orbital t a,b show inverse bifurcation at the 
FM-FMz phase transition being the end of the AO order but remain non-zero in the paramagnet 
which is induced by finite E z . The on-site SOE represented by r a ' b are small everywhere but in 
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Figure 2. Panel (a)- phase diagram of the d 9 monolayer at T = 0; for a more complete phase 
diagram see [13]. Shaded (yellow) area marks singlet phase with spin disorder. Panel (b)— spin 
s and orbital t a ' b order parameters, on-site r a,b and bond R a,h SOE (multiplied by 100) for 
point (E z = —2.1J,7] = 0.19), marked with red arrow in panel (a), as functions of temperature. 
Thermal evolution of the system involves intermediate FMz phase before reaching PM phase. 



FMz they become relatively large which lets us think of FMz as an entangled phase. The bond 
SOE defined by R a,b = (SiSjT""' r"' 6 ) — (SiSj)(r] 2 ' 6 r"' b ) with (li)\\a(b) take considerable values 
only in the PM phase and shows rather purely statistical covariance than a true quantum effect. 

Summarizing, we have shown that spin-orbital entanglement leads to exotic types of order 
which are stabilized by quantum fluctuations in a bilayer and by thermal fluctuations in a 
monolayer. They emerge from highly frustrated spin-orbital superexchange and could be 
discovered only within a cluster mean field approach. 
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